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Abstract 
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the weak solutions to the primitive system converge to the solution to the 
anelastic Navier-Stokes system and the transport equation for the second order 
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1 Introduction 

We study certain singular limits arising in the analysis of stratihed fluid flows in 
meteorology and astrophysics, mm- Moving beyond the case of homentropic 
flow, which was studied rigorously by Masmoudi, [19], we include variations of an 
entropy variable - the potential temperature, 0, at second order in the Mach number. 
This is motivated by the observation that only when such entropy variations are 
accounted for, will the fluid flow equations support internal gravity waves, and these 
are responsible for a host of important physical processes, see, e.g., [T]. Specihcally, 
we study the asymptotic regime introduced originally in a formal asymptotic analysis 
by Ogura and Phillips [22]. Within this regime, advection and internal waves act on 
comparable time scales while sound waves are asymptotically fast in the low Mach 
number limit. Thus, considering this flow regime allows us to incorporate the effects 
of entropy stratihcation without the added difficulty of an asymptotic three-scale 
problem, [15], that would have advection, internal waves, and sound act on three 
asymptotically separated time scales. 
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An additional simplifying assumption concerns the time evolution of the poten¬ 
tial temperature. In the mentioned application areas, some mean potential tem¬ 
perature stratihcation is generally observed over long time scales which differs from 
that associated with constant temperature, T. These stratihcations are maintained 
against dissipative processes that tend to homogenize the temperature by various 
mechanisms including radiation, turbulent motions associated with convective in¬ 
stabilities, latent heat conversion, and similar “diabatic effects”. In the present hrst 
approach to addressing flows with potential temperature stratihcation, we exclude 
the mathematical complications induced by these processes and instead ensure the 
maintainance of a mean potential temperature stratihcation by excluding diabatic 
processes and the inhuence of molecular transport on potential temperature, and by 
imposing a mean stratihcation in the initial data. 

Under these premises, the huid motion is described by the Navier-Stokes system 
in the isentropic regime. Specihcally we analyze the following dimensionless system 


of equations for e 1: 

(9tp-Fdiv3,(pw) = 0, (1.1) 

dt{gu) + diva;(pn 0 n) -f ^V 3 ;(p 0 )'^ = udiv^E>(yxU) + (1.2) 

dt{gQ) + diva,(p 0 w) = 0 , ( 1 . 3 ) 

where g denotes the mass density, u the velocity, and S(V 3 ;w) the viscous stress 
tensor, here given by Newton’s rheological law 

S(Va;'a) = p (VxU + — -diva;itlj -|- Adivajul, p > 0 , A > 0 . ( 1 . 4 ) 


Here we have adopted a distinguished limit that equates the Mach and (external 
wave) Fronde numbers with the small parameter, i.e., Ma = Fr = e, while the 
Reynolds number Re = 1/u is treated as an independent parameter. In part of 
this work, we consider limit processes, in which e —)■ 0 and z/ —)■ 0 independently. 
The Strouhal number is unity as we adopt the advective time scale tref = ^ef/wref 
as the reference time scale for non-dimensionalization and, following the discussion 
of the last paragraph, we set the Prandtl number, Pr = cxo, thereby neglecting the 
molecular transport of heat. 

The fluid occupies a slab 

n = nhX{o,i), (1.5) 

where the horizontal projection Qh is either a periodic “flat” torus T 2 D = ([0, l]|{o,i})^ 
or the unbounded space M^. The velocity is supposed to satisfy the complete slip 
closed lid boundary conditions. 


u n = uslon = 0, (S(V^w) ■ n) x n|aQ = 0, (1.6) 

where n denotes the outer normal vector to dQ. We shall denote x = [xh, z], where 
Xh = [xi,X2] 

We take F = —gz, g > 0; whence = [0, 0, —g] represents the effect of the 
gravitational force acting on the fluid in the vertical direction. Accordingly, and 


3 


since we assume the potential temperature to be constant to leading order (see fll.Sp 
below), the leading order equilibrium distribution of the density g satishes 


= Q^xF. 


(1.7) 


Obviously, equation fll.7p admits a continuum of solutions, here we suppose g = g{z) 
to be given, g{z) > 0 for all 2 ; G [0,1]. 

Problem fll.lHl.6l) is supplemented with the initial data 


^(0,-) = ^o,£ = ^ + 

< ^( 0 , •) = > 


( 1 . 8 ) 


. 0(0, e) — 0o,£ — 1 + , 

Moreover, if VLh = is unbounded, the far-held conditions 

g ^ g, 0 —)■ 1 as |a;| —)■ cxo 


(1.9) 


are prescribed. 

Our goal is to study the behavior of solutions [g^, u^, O^] in the asymptotic limit 
£ —)■ 0. Formally, it is not difficult to check that 

~ ^ as e —)■ 0 , ( 1 - 10 ) 

whereby fll.ip reduces to the so-called anelastic constraint 

div3,(^n) = 0 (1.11) 

for the limit velocity v. Moreover, the specihc form of the initial data for 0 indicates 
that 

0 ( 2 ) ^ ^ ^ T for £ —)■ 0 , ( 1 - 12 ) 

where T is transported by v, 

atT+n-V^r = 0. (1.13) 

Finally, we also perform the limit in the momentum equation fll.2p . Assuming that 

also the coefficient z/ —)■ 0 vanishes in the asymptotic limit we recover, again formally, 
the system 

dtv + V-VxV+ Va:Il =-TVa:F, (1.14) 

where fl is the pressure, the presence of which is enforced by (11.111) . In the case 
z/ > 0 is kept hxed, a viscous tensor |§(Va;'n) will necessarily appear in the limit 
system. 

We consider weak solutions to the primitive system fll.lH1.3l) according to the 
dehnition given in section [2l The advantage of such an approach is that these 
solutions are known to exist globally in time, the drawback is their low regularity 
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given by the typically very poor a priori bounds. Two kinds of initial data will 
be considered: (i) well-prepared data mimicking the structure of the data of the 
expected limit system, (ii) ill-prepared data which merely require uniform bounds 
for certain rescaled quantities as e —)■ 0. We show that in both cases the weak 
solutions [pg, u^, ©g] of the primitive system converge to the expected limit. In the 
case of the well-prepared data, we also obtain an explicit rate of decay in terms of £ 
and u assumed to vanish in the asymptotic limit. The method used is that of relative 
energy (entropy) providing a suitable “distance” between the weak solutions of the 
primitive system and that of the target system. Accordingly, the convergence takes 
place on the life-span of (smooth) solutions to the limit problem. Ill-prepared initial 
data give rise to high-frequency acoustic waves that must be hltered out in the limit. 
To this end, we consider the unbounded slab hi = x (0,1) allowing for dispersion 
of acoustic waves that actually vanish in the asymptotic limit on any compact subset 
of n. The analysis at this point relies on the application of the celebrated RAGE 
theorem. In contrast with the case of well-prepared data, we establish global-in¬ 
time convergence towards a weak solution of the limit problem. Unfortunately, as 
uniqueness of the weak solutions for the limit problem is not known, there is no 
explicit rate of decay available for the ill-prepared data. 

The paper is organized as follows. In Section [2l we recall the concept of weak 
solution for both the primitive and the target system used in this paper and state the 
main results. In Section 121 we introduce the relative (modulated) energy functional 
needed in this context and derive the necessary uniform bounds independent of 
the scaling parameter e. Section 0] is devoted to the case of well-prepared initial 
data while Section |5] deals with ill-prepared data. Possible extensions and further 
comments on the methods used in the paper are discussed in Section 01 

2 Weak solutions, main results 

In this section, we introduce the concept of weak solution to the Navier-Stokes 
system fll.lHl.Sp and state our main result. 

2.1 Weak solutions to the primitive system 

We consider the weak solutions to system fll.lHl.Sp belonging to the class 

g E GweakdO,T]; T"^(hl)), 

0 e L°°((0,T) X f]), (2.1) 

u E 

and enjoying further regularity and integrability properties allowed by the available 
a priori bounds. We remark that the impermeability constraint 

= 0 ( 2 . 2 ) 

makes sense in the above specihed class. 
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2.1.1 Weak formulation of the transport equations 

We say that u is a weak solution to the equation of continuity fll.ip in (0, T) x Q 
if the integral identity 


gip da: 


t=T2 


t=Tl 


[gdtif + gu ■ Vx^] dx dt 


' T\ J Q 


(2.3) 


holds for any 0 <Ti <T 2 <T and any test function (p G C“([0, T] x 12). 
As for equation fll.dp . we consider its renormalized version, specihcally, 


gG{Q)ip dx 


t=T2 


t=Tl 


[gG{Q)dtp> + gG{Q)u ■ Vx'p] dx dt 


(2.4) 




for any (p G C“([0,T] x 12) and any G G C(M). 

Remark 2.1 Alternatively, we can replace ^2.4\ ) by postulating the renormalized 
version of the pure transport equation 


dtO + u ■ Vx0 = 0, 


(2.5) 


namely 


G{Q)ip dx 


t=T2 


^ J t=Ti ari JQ 

for any ip G C))°([0,T] x 12) and any G G C'(M) 


[G{Q)dt(p + G{Q)u ■ VxP^ + G(0)div3;U(/?] dx dt 

(2.6) 


2.1.2 Weak formulation of the momentum balance, energy inequality 

The standard weak formulation of the momentum balance fll.2p reads: 

r 1 t=T 2 


gu ■ ip dx 


Jei 


t=Tl 



gu ■ dfip + gu^u : VxP> + -^(p 0 )^div 3 ;(p 


ri O 


(2.7) 




dt 


for any ip G C“([0, T] x 12; M^), ip ■ n\QQ = 0. 

In addition, we focus on the class of finite energy weak solutions satisfying the 
energy inequality 




1 , ,2 1 

+?(7VT) 


< 


(p0)' 


dx 


t—T pr p 

+ 1x1 / 'Bi{VxU) : VxU dx dt 

t=0 Jo 


(2.8) 



;gVxF ■ u dx dt. 


0 Jn 
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It was shown by Michalek [21] that the family of weak solutions to problem fl2.3l) . 
fl2.6p . fl2.7p . and fl2.8p is weakly sequentially compact, meaning any sequence of weak 
solutions contains a subsequences weakly converging to another weak solution to 
the same problem. Novotny et al. mi established existence of global-in-time weak 
solutions in a slightly different setting for the unknowns [p, u, Z = p0], where f|2.4p 
is replaced by 

dtZ + (\iYx{Zu) = 0. (2.9) 

For further comments concerning the existence of weak solutions satisfying fl2.4p 
and/or fl2.6p see Section [6l 


2.2 Solutions to the target system 

As observed in the introductory part, the target system in the inviscid limit is fll.llF 
I1.14p . In view of the global existence result by Oliver [23l Theorem 3] for 2D system, 
we may anticipate the existence of local-in-time strong solutions [n, T, 11] also in the 
3D case, specihcally. 


n,r e C'([0,T^ax);hF"*’'(f^)) 


(m > 3) 


( 2 , 10 ) 


provided 


t;(0,.) = no e 

r(0,-) = roelF™’ 2 (D), 

^^0 ■ nlsn = 0, 


( 2 . 11 ) 


diVa;(pno) = 0. 

For hxed positive viscosity coefficient u > 0, the limit problem is the Navier-Stokes 
system 


div,^(pn) = 0, 

gdtv gv ■ + pV^Il = i^div3,S(Vx'y) - gTVxF, 


( 2 . 12 ) 

(2.13) 


where T satishes the transport equation (I1.13p . The weak solutions of (I2.12p . (I2.13p . 
supplemented with the slip boundary condition (II.6p are dehned in a standard way 
by requiring the integral identity 



t=T2 

t=Tl 


(2.14) 



gv ■ dtip + [gv Z>v - uE>{Vxv)^ 


Vx(p + gTVxF ■ if 


dx dt 


to holds for any (p G C/°([0, T] x D; M^), diva;(pv9) = 0, (p ■ n\QQ = 0. 
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2.3 Main results 


Having collected all the necessary preliminary material we are ready to state the 
two main resnlts of this paper. We start with the inviscid limit for the well-prepared 
initial data. 


Theorem 2.1 Let 12 = T 2 D x (0,1) and 7 > §■ Let ©e] he a weak solution of 

the primitive system U.l[\1.6\) on a time interval (0,T), with the initial data U.8\) . 
where 

pg, ©g e L-(12), e L“(12;M=^), (2.15) 

ll 2 oilU°°(o) + l|0oilU°°(f2) + II'^0,£|1 l°°(O;R3) < D. (2-16) 

Suppose that the target system U.lljfTJff j admits a smooth solution [r;, T] on the 
same time interval [0,T] emanating from the initial data 

v{0,-) = vo, T{0,-) = To. (2.17) 


Then 


snp 

te[o,r] 



Qslu^ - v\'^ + 



7 

+ Qe 



dx 


( 2 . 18 ) 


<c(T,D) 


£ + V+ / |'Uo,e~'*^o| + 


2 o,£ 


+ 


©( 2 ) 


Tn 


dx 


I Jn 

where the constant c = c{T, D) depends on the norm of the limit solution [v, T] and 
on the size D of the initial data perturbation. 


Remark 2.2 Theorem \2.1\ yields convergence in the regime e —)■ 0, z/ —)■ 0 and 
for the well prepared initial data for which the expression on the right-hand side of 
l[2.1S\) vanishes in the asymptotic limit £, o 0. 


©nr second resnlt concerns the case of ill-prepared data. 

Theorem 2.2 Let 12 = x (0,1) be an infinite slab and 7 > 3. Let [g^, ©g] be 

a weak solution of the primitive system on a time interval (0,T), 

with the initial data M.^) . where 

So 5 

To weakly-(*) m L^ n L°°(12), (2.19) 

Vo weakly-(*) in Lf fl L°°( 12 ;R^). 

Then, passing to a suitable subseguence as the case may be, we have, for e 00 , 

SUPig(o,T) Ike - 2|U2 (q) ^ 0, 

< Us ^ V weakly in L'^{0,T-,W^’‘^{Ll-,M.^)) and (strongly) in Ly^^{[0,T] x > 

7 " weakly-(*) in L°°{{0,T) x 12), 


( 1 ) 


2o, 


-)■ 


©( 2 ) _v 

U0,e 


( 2 . 20 ) 



































where [i>, T] is a weak solution to the target system H2.1B) . ([ 2 . 14 ^ , with the 

initial data [i;o, To]- 


The rest of the paper is basically devoted to the proofs of Theorems 12.1112.21 


3 Relative energy, uniform bounds 


Similarly to Masmoudi [18], Jiang and Wang [2l], among many others, we adapt the 
idea of Dafermos [1| based on the concept of relative (modulated) energy associated 
to the compressible Navier-Stokes system. In the present context, the relative energy 
functional reads 




0,-u r, it) = f ^g\u-U\'^ + - H\r){ge-r) - H{r)^ 


where we have set 


H{Z) = 


7 - 1 


dx, 

(3.1) 

(3.2) 


3.1 Relative energy inequality 

As shown in 1131. any finite energy weak solution to the Navier-Stokes system satishes 
the relative energy inequality 




r,U 


t=T 


J t=0 



+ u / §(V^(w-C/)) : V^(w-C/) dxdt 


0 n 



< / g{dtU + u-VM)-{U-u) + nS{VM)-.^AU-u) 


0 o 


(3.3) 


+-{{r- ge)dtH\r) + V.JJ'(r) • {rU - gOu) 


- divM ({gOy -r^) - gV^F ■ {U - u) } dx dt 


for all (smooth) “test functions” functions r, U, 

{r-g)e ^“([0, T] X H), U e ^,""([0, T] x H; R=^) (3.4) 

such that 

r > 0, U ■ n|an = 0. (3.5) 

Remark 3.1 A special form of liS. A) for a particular choice of the test functions r, 
U and under different assumptions on smoothness of the solution [p, u, 0] was de¬ 
rived by several authors, e.g., Desjardins Germain Mellet and Vasseur fWf. 
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3.2 Uniform estimates based on the (relative) energy in¬ 
equality 

The choice r = g, U = 0 in the relative energy inequality fl3.3p gives rise to 


^Qe\ue\‘^ + - H\g){g^e^ - q) - H{g)^{T, •) 

S{VxUe) ; VxUe dx dt 


dx 


< 


in . 


+v 

1 



0 Jn 


^*0,e|’^0,e| T ^2 H (^) (^0,6:00,e H 


dx 



0 Jn 


VxH'{g) ■ dx dt + — / / gsV^F ■ dx dt. 


0 Jn 


In view of fl3.2p and 


VxH\~g) = VxF 


(3.6) 


we may infer 


In . 


< 


in 


+v 

1 



0 Jn 


< 


Qe\ue\'^ + -^{H{geQe) “ H\g){gSe -q)- hf(^))(r, •) 

E)(VxUs) : dx dt 

2 ^ 0 ,e|'t^ 0 ,eP + ^ {Q 0 ,e^ 0 ,e) ~ (^ 0 , 6 : 00 ,e — Q) ~ H [g) 

\ [ [ ge'^xF ■ U^{1 - Qs) dxdt 
^ Jo Jn 

^*0,£|”^0,e| T ^ ^.^(^*O,£0O,e) H (^)(^O,£0O,£ 


dx (3.7) 


dx 


dx 


IV.FII 


L°°( 0 ) 



0 Jn 


Qe \ 'U'e \ T Qe 


0 . - 1 , , 

dx 


dt. 


Finally, taking 9 ? = 1, G = ( ^^2 ^ ^ in (12.4p we obtain 


Qe 


0 £-l 


) (t, ■) da: < J go,, (©S) 


(3,8) 


Relations (13.7p . (13. 8 p yield uniform bounds on the family {^g, Q^, ite}£>o provided 
the initial data are bounded as in the hypotheses of Theorems 12.1112.21 specihcally, 
are chosen in such a way that 


|'n'0,£||LinL°°(O;R3) < C, 


(3.9) 
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Qo,e = Q + lko^el|LinL°°(0) < C, (3.10) 

0o,£ = 1 + ^^0oi) ||0oilUinL°°(Q) ^ c. (3-11) 

Now, we deduce from fl3.7l) . (13.81) the following uniform bounds independent of e, 


V. 


iG[0,T] 


VxUe + - -diV2;tt£l 


ess sup 

iG(0,T) 




ess sup / 

tG(0,T) Jg/2<geee<2g 


i?£ 0 £ Q 


L^{n) 
2 


C, 

(3.12) 

< c, 

,r)xO) 

(3.13) 

< c, 

1 

(3.14) 

dx < c. 

(3.15) 


ess sup / 

tG(0,r) Jg,ee<e/2 


(1 + |p£ 0 eP) dx + ess sup 


te(o,T) JQ^e^>2e 


(1 + |p£0£r) dx < e^c. 


(3.16) 


As a matter of fact, relation ( 12 .4p may be used repeatedly to deduce 

< c (3.17) 


0 £-l 


as well as 


ess sup 

46(0,T) 


ess sup / Og 
4G(0,T) Jn 




^ dx < c. 


(3.18) 


Remark 3.2 As the function appears in only multiplied by g^, we 

may assume, without loss of generality, that ©^ = 1 on the (hypothetical) vacuum 
set pe = 0. Keeping this convention in mind, we observe that (jS-lTO holds in the 
whole physical space hi. 


Now, we have 


[ 

Qs P 

Jg/2<geee<2g 

e 


dx < 


/ 

Pe(©£ - 1) 

dx+ [ 

P£0£ P 

) gl2<geQe<2Q 

£ 

Jg/2<geee<2g 

e 


dx. 


therefore (I3.14p . (I3.15p give rise to 


ess sup 

4G(0,7 

Similarly, the bounds (I3.16jl . (I3.17p yield 


f 

Qe P 

}g/2<g^Be<2g 

e 


dx < c. 


ess sup 


|Pepdx + ess sup 


I Pep dx < e^c. 


te{0,T) Jg,ee<g/2 


t€(0,T) Jg,ee>2g 


(3.19) 


(3.20) 
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At this stage it is convenient to introdnce the following notation; for each 
measurable function f(t,x) we write / = [/]ess + [/]res, where 


[/]ess = [/]res = (1 “ 


(3.21) 


X = 


1 for g/2 < Qe^e < 2^ 
0 otherwise. 


Now, we have 


e,-i 


da; < 


Q - 




0e-l 


da; + 


Qs^e 


0e-l 


da;; (3.22) 


where the last integral is bounded due to (13.17p and fl3.18p . As for the first integral, 
we may decompose the integrant into its essential and residual part, where the 
essential part can be estimated as follows 


2e0e 


0.-1 




Q 


0.-1 


0.-1 


da: 


da:. 


while for the residual part we have 


(3.23) 


/o L 


2.0. 


0.-1 1 

2 

0.-: 


£2 

f 

£2 

L°°{n) J o 


da: 


Q - 


2 e 0 e 


da; < c 


(3.24) 


which follows from estimates (I3.16p and (I3.17p and the fact that ||p||l°o(o) < c. 
Putting together estimates (13.22113.241) and recalling that p > 0, we deduce 


ess sup 

*6(0,T) 


0.-1 


< c. 


(3.25) 




4 The inviscid limit 


In this section, we prove Theorem 12.11 To this end, we take C/ = v, r = p in as test 
functions in the relative energy inequality (13.3p . where v is the (smooth) solution 
of the target system (ll.lUfTTTp . After a straightforward manipulation (see [11] for 
more details) we obtain 




^ 2 .) 0 .) 


u. 


g,v 


t=T y 

J t=Q 2 


' Tl J Q 


E>(y^u^) - s(y^v) 


WxUe — 'VxV ) da; df 


(4.1) 



< / / 2. {dfV + ■ Vxv) ■ {v — u^) dx dt + v I / | da; dt 

'0 Jn Jo Jn 
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0 -'O 



- QeQe)dtH'{Q) + Va;H'{g) ■ {gv - g^Qe^ 

1 


0 Jn 


div^'y(^(^e0£)^ - dx dt- — 



dx dt 

Qs'^xF ■ {v - u^) dx dt 


0 Jn 



< ||Vxi^||L°°((0,T)xn;R3x3) / / ge\us-v\ dx + u i / \Vxv\ dx dt 


0 Jn 



0 Jn 



>0 Jn 
1 



ge {-TVxF - V^n) ■ {v - Ue) dx dt 

'^xH'ig) ■ (gv - gsQeUe) dx dt 
1 



0 Jn 


0 Jn 


- {gTy^ dx dt - — 




0 Jn 


II II L°°^^g j'^xQ.jjSxS^ / / 1^1 dx-\- ly 


Qe^xF ■ {v - uy da; dt 
IVa^t^l^ dx dt 


0 Jn 




0 Jn 


Qe ( -T-^ VxF ■ {v - uy dx dt + 



0 Jn 

• (n^ - v) dx dt 


0 Jn 


+— / / VxH\g) ■ {gv - gsQeUe) dx dt - - / diva;a;('(^e0£)'^ - dx dt. 

£ ./o Jn ^ Jn Jn X / 


0 Jn 


Next, we use (I3.6p to rewrite the penultimate term as follows 

'^xH'{g) ■ {gv - geOeuy dx dt 
1 



0 Jn 



0 Jn 



{g - geQE)'^xH'{g) • u dx dt + ^ / / geQe^xH'{g) ■ {v - Ue) dx dt 


0 Jn 



{g - geQe)'^xH'{g) ■vdxdt + 


0 Jn 



gE^e^xF ■ {v - Ue) dx dt. 


0 Jn 


Furthermore, in accordance with the anelastic constraint fll.lll) . 

Vxg-v = -gdivxv; 

whence 

{g - gE^y'^xH^g) ■ v dx dt = \ [ [ (^£0^ - ^) 7 ^'^"MiVa,t; dx dt. 

0 Jn ^ 



0 Jn 


Thus, revisiting fld.ip . we deduce 


g,v 


t=T 


0£) 

< l|V.i;|| L°°((o,T)xr2;R3x3) 



§(Va;We) - S(Va:'u)) : ( 


VxUe — Vx'U ) dx dt 
(4.2) 


uj dxdt + 12 / / dx dt 

J Jq Jn 
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06 


0 Jn 


0,-1 


— T I 'VxF ■ (v — Ui;) da; dt + 



■ (u, - v) dx dt. 


0 Jn 


The last integral on the right-hand side can be handled as follows: 

n PeVa;!! • (w, - v) dx dt 

I 



Pedfll dx + 


0 ao 


^,11 da: 


t=T 


J t=o -'0 Jn 


+ {0- ■ V da:; 


whence 



• (-u, - v) da: dt 


'0 Jn 

Going back to fl4.2p we obtain 


< c ess sup Up, - p||li(o)- 
te[o,T] 


(4.3) 


^6 \ 061 0e) 


01 V 


t=T 


u=o 


< C (II V3;r^||L°°((0,T)xO;R3x3) + || 2,oo((o^r) xO;R3)) 

iVri^P da: dt 



0 JQ 


^6 [061 0e) '^e 


(4.4) 
Q, dxdt 


+1^ 



0 Jn 


+c||V,F|U °°((0,T)xO;R3) ' / 06 



0 Jn 


0,-1 


T ] dx dt 


-fc||n||^i,oo((o^r)xO) ess sup ||Pe — pIIlPo). 

tG[0,T] 

Finally, to bound the penultimate term we show: 

Lemma 4.1 Suppose g, Q is a weak solution of ft.,91) in the sense specified 

dtT+v- VxT =0, T eC\ve C\ 


in ^2.4\). Let 


(4.5) 


Then 


-p|G(0)-rp da: 


Proof: 

We write 


t—r 2 pT 2 r 

- t=Tl 


g(^G{Q) — (^u — ■ WfiT da: dt. (4.6) 


-p|G( 0 ) -n^ da: 


1 


pG(0)2 dx 


Lao 

t=T2 


+ 


t=Ti 


gG{e)r dx 


t=T2 

t=Tl 

t=T2 


J t=Tl 


-gT^ dx 


t=T2 


J t=Ti 


where, in accordance with fl2.4l) . 


Lao 


]^gG{Qf dx 


t=T2 


= 0 . 


t=Tl 
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Next, using the weak formulation fl2.4p with the test functions (p = T and again 
with G = 1, if = we obtain 


gG{e)T dx 




t=T2 


-qT^ da; 


t=T2 


t=Tl 


t=T\ L'J 

QG{Q)dtT + qG{Q)u ■ V,r + gTdtT + gu • TV^T 


J Tl J Q, ^ 

whence, by virtue of fld.Sp . 


da; df; 


gG{Q)T da; 


.Jn 


t=T2 


t=Tl 


-gT^ da: 


In 


t=T2 


t=Tl 


^ Tl J Q. ^ 


— gG{Q)v ■ VxT + gG{Q)u ■ VxT — gTv ■ VxT + gu ■ TVxT 
= f" f 0 ( 0 ( 0 ) -r](u-v] ■ VxT da; dt. 


da; dt 


^ Tl J Q, 


Q.E.D. 

Using (14. 6 p and the uniform bounds established in (13.124134(711 . we may now 
apply Gronwall’s lemma to (I4.4p . to obtain convergence towards the limit system. 
In addition, the rate of convergence can be computed explicitly in terms of e, v as 
indicated in fl2.18p . Theorem 12.11 has been proved. 

5 The ill-prepared initial data 

With u > 0 hxed, our ultimate goal is to preform the singular limit e —)■ 0 for 
the ill-prepared data as claimed in Theorem 12.21 This is dehnitely more delicate 
than the previous part as the motion is polluted by rapidly oscillating acoustic 
waves. Fortunately, the unbounded physical spaces gives rise to dispersive effects, 
in particular, the acoustic part of the velocity vanishes in the asymptotic limit, at 
least on compact subsets of 42. 

5.1 Weak convergence 

To begin, a generalized version of Korn’s inequality (see e.g. m Section 10.9]) can 
be used to deduce from fl3.13p . fl3.19p . and fl3.20p the bound 

I \VxUe\‘^ + \u^\‘^ dx < c, (5.1) 

n 

in particular, we may assume that 

^ V weakly in L^(0, T; W^’‘^(VL] M^)), 
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passing to a subsequence if necessary. Similarly, it is easy to check that 

ess sup II Pe - p||l2(0;R3) 0 

*6(0,T) 

keeping in mind that 7 > 3 in Theorem 12.21 

Finally, it is easy to observe that the continuity equation fll.ip reduces to the 
anelastic constraint fll.lip in the asymptotic limit e ^ 0. 


5.2 Momentum balance 

We write the momentum equation 


in the form 


dtiQeUe) + - - Pe0.V.F) 

^ 10 

= ® Us) + udWr,S{VxUs) + Qe - 


(5.2) 


where, furthermore, 

- Qe^e^.F 

= [(PeOs)'^ - -Q)- Q^] 

-lQe^eQ^~‘^'^xQ + - ^)] 

= Vx [(peOe)'^ - lQ^~^{Qe’de - q) - f]-{QeQe-Q)lQ^~‘^'^xQ+'^x [lQ^~^{QeQe “ ^)] 
= [(peOe)'^ - lQ^~^{,Qe^e - Q) - f] + Q'^x {Qs^s “ ^)] • 

Consequently, a suitable weak formulation of (15.21) reads 

/ / [qsU^ ■ dtip + QeUe ® : Wx^] dx df 


(5.3) 


Jo Jn 

f {QeQsV - 10^~^{QeQe - Q) - j 

- aWxV dx df 


'0 JQ. 




'0 Jn 


f \ 1 — 0 

/ / V^{VxUe) : Vx^ - Qe - xF ■ v? 

'0 io L ^ 


— ^ [qSs - q)) div3,(pv9) dx df 

dx df - / Qo,eUo,e ' <^(0, •) dx 


for any ip G C“([0,T) x Vt, ip ■ n\gQ, = 0. 

Next, in agreement with the uniform bounds (I3.17p . (I3.25p . we may assume that 


0.-1 


^ r weakly-(*) in L°°(0, T; n L“(f2)). 


Consequently, it is possible to perform the limit for e —(■ 0 in both (15.3p . for ip 
satisfying div 3 ;(p(p) = 0, and (II.3p to obtain 


'0 Jn 


[gv ■ dtip + gu ® u : Vx^] dx dt 


(5.4) 
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lim r f 


€—^0 


'0 Jn 


= / [z/§(VccV) : - qTVxF • v?] dx dt- / qvq ■ <^(0, •) dx 

Jo Jn Jn 

for any G C“([0,T) x f2, (^ ■ n\gQ = 0, div 3 ;(^(^) = 0, and 


/ / [QTdt^p + qTv ■ VxV?] dx dt = - / ^ 7 oV 5 ( 0 , •) dx 

'0 io in 


(5.5) 


for any ip G C‘^{\^,T) x f2), where qu®u denotes a weak limit of the sequence 
{Q^Ue®Ue}e>0- 

Remark 5.1 Note that we have used Lions-Auhin lemma to obtain / (5.5)) . 

Thus, in order to complete the proof of Theorem 12.21 it remains to show: 


lim r f ^ (5 g, 


£—^0 


0 JQ, 


QUt^U=QV(SV. 


(6.7) 


This will be done in the next section by means of careful analysis of propagation of 
acoustic waves. 


5.3 Acoustic equation 

We start by introducing a generalized Helmholtz decomposition 

w = Vg{w) + gQg{w), Qg{w) = Vxd/, (5.8) 

where T is the unique solution of the Neumann problem 

div 3 ;(pV 3 ;T) = div^-iu, {w — ■ n\gfi = 0 , T —)■ 0 as |x| ^ oo. 

Accordingly, we write 

QeU^ = VgigeUe) + gQg{geU^)- 

Note that, by virtue of (13.121) . (13.201) . (15.ip . and the hypothesis 7 > 3, we may write 

ge'^e [^?e”^£]ess T [^?e'^£]res) 

cf. (I3.2ip . We have 

[^£H<.]ess -t gv weakly-(*) in L°°(0, T; L^(H; M^)), 
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while 

Seeing that 


[Qe^elres “t 0 in L^((0,T) X 


w ■ if dx = / Vq{w) ■ ip dx whenever div 3 ;(^tn) = 0, 


we may use the momentum equation fl5.3p . together with the bounds established in 
Section [221 115.ip . and the standard Lions-Aubin argument to conclude that 


'Pg{QeU^) VgiQv) = QV (strougly) in Lf„^((0,T) x 
To derive an equation for the acoustic potential we introduce notation 

Qe^e Q 


(5.9) 


S,= 


eg 


— Qg^Qs'^e') • 


First, we rewrite equation fll.3l) in the form 


sdtS^ + -div3,(pVj;<F£) = £-div,j, g^ 

g g \ \ e 


1 - 0 , 


u. 


whereas 05.31) reads 


= -eV, 


edtigsu^) + pVa; [yg'^ 
- 7Q'^~\QeQe - Q) - 


(5.10) 

(5.11) 


ediVooigeUe ( 8 ) Ue 


+si^dw^^{V^u^) + age -— 


Our goal is to apply the generalized Helmholtz projection to 05.lip . This 
step may be performed in the weak sense by taking a test function 

e C“((0, T) X H; M^), ■ nlan = 0 

in the weak formulation of 05.lip : 


'0 Jn 


[agV^^e ■ dtV^cp + div3.(pV,j0)] dx dt (5.12) 


= —£ 


'0 Jn 


{QeQeV - 10^ ^{Qe^e - Q) - 0^ 


A0 dx df 


—£ 


72 J. ,,(S37V7 ^ . V72 J. , „ 1 —V F ■ V 


'0 Jn 


{geU, (g) u,) : - I/S(V^'U£) : V^(/) + g,- 




dx df. 


Now, for ip G C“((0, T) x H; M^, ip ■ nl^Q = 0, we may take Vx((> = Qgiyj) as a test 
function in 05.12p to obtain: 


'0 Jn 


[-e^edtddvx,ip -f (yp'^ div^-y:] dx dt 


(5.13) 
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= —e 


/o Jn L 


- q)- 


diVj.Q^((^) dx dt 


—e 


'0 Jn 


[{QeUe ® Ue) : VxQgi^) “ V^{Vx^e) : VxQq{^)\ dx dt 


—e 


0£- — 5 ^V^F • Qgi^p) dx dt 


Jo Jn 

We conclude that the system (15.101) . (15.131) can be written as a variant of 
Lighthill’s acoustic analogy: 


edtS^ + -div^(^V,:<h£) = eQ], 

Q 

edtVx^e + ^. {rQ^-^S,)=egl 
supplemented with the Neumann boundary condition 


• 'M'lgo — 0 , 


(5.14) 

(5.15) 

(5.16) 


where we have set 


<3l = -div^ ( pe 


1 - 0 , 


u, 


and 

Ge = “QgVx 


- Q)- 


1 - 0 , 


-Qpdiv^(p£'U,®it,)-Qgdiv^S(V^ite 


+£Q50e - 5 - '^xF. 


5.4 Analysis of the acoustic system 

Consider the homogeneous problem associated to (15.141) . (I5.15p . namely 

dtS + 4div^(pVx$) = 0, dt^ + c(^)S' = 0, Vx^ ■ n\on = 0, (5.17) 

Q 

where we have denote 

c(p) =p'{,q) = 

that can be viewed as a wave equation 

52 ^$ _ £i|ldiv^(pV^<l)) = 0, V^d) ■ n\QQ = 0 
0 

for the potential < 1 >. 

The acoustic propagator 

A{w) = -c{ 0 )AyW - { 0 dzw), X = {yi, 2 / 2 , z), VxW ■ n\on = 0, 

Q 
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has been studied by several authors. In particular, DeBievre and Pravica 0 . 0 
showed that A can be viewed as a self-adjoint operator on a weighted space 
endowed with the scalar product 


(u,v}^ = 



Q 



dx. 


Following the strategy of [9] we derive the desirable dispersive estimates for the 
solutions of the acoustic system fl5.14l) . flS.lSp from the following version of the 
celebrated RAGE theorem (see Cycon et ah [3], Theorem 5.8]): 


Theorem 5.1 Let H he a Hilbert space, A : ^[A) <Z H ^ H a self-adjoint operator, 
C \ H ^ H a compact operator, and Pc the orthogonal projection onto the space of 
continuity He of A, specifically. 


H = Hr 


clj:/|span{tc E H \ w an eigenvector o/A}|. 


Then 


T 


exp(—itA)CPc exp(hy4) dt 


—0 as r —oo. 


(5,18) 


cm 


5.4.1 Absence of the point spectrum for the wave propagator 

In order to apply Theorem 15. II we have to make sure that the point spectrum of A is 
empty. To see this we use the positive commutator method introduced by DeBievre 
and Pravica [B] in a similar context. More precisely, we compute the commutator 

[A, y ■ Vy] (w) = A{y ■ Vyw) - y ■ Vy(A(M;)) (5.19) 

= -c(p) [Ay{y ■ Vyw) - y ■ Vy{Ayw)] = -2c{g)AyW. 

If 

A{w) = Xw, 

then, since A is self-adjoint, we get 


/ w[A,y ■ Vy] w dy dz = / {A{w)y ■ VyW — Xwy ■ '^yw) dy dz = 0, 

Jn Jn 

while, in accordance with 05.191) . 


w[A,y ■ Vy] w dy dz = 2 / c{g)\'VyW 


J fi Jfi 

Thus any possible eigenfunction w must be constant with respect to the horizontal 
component; whence ta = 0. 
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5.4.2 Decay of acoustic waves 

Following step by step [9l Section 5] we rewrite the acoustic system fl5.14p . fl5.15p in 
the abstract form 

edt{Z„(p)^+ = (5.20) 

for any smooth 0 G Va;0 • 'M|ao = 0, 

edt ($£, <P)h + {Ze, <P)h = ^ ^)h ’ smooth 0 G (5.21) 

where we have set = c{g)Ss. The quantities are computed in terms of 

Q^. After a bit tedious but straightforward computation for which we refer to 
[9l Section 5], it can be shown that 

7l = H\A)\h% Tl = E\A)\h% 

where 

{hl}e>o, {hi }e>o are bounded in L^((0,T) x f2), 

and ; (0, cxo) M are suitable smooth functions that may become singular 

at 0 and cxo. 

Now, exactly as in [9l Section 5.2], we apply Theorem 15.11 in the situation 
A = VA, C = x"G{A), X G cr (H), X > 0, G G Gr (0, oo), G > 0, 


to deduce that 

II<^(-4)[^£]|Il2((o,t)xA') -t 0 as e 0, (5.22) 

||G(.A)[Z£]||^ 2 ((o,T)xi^) -t 0 as e -)■ 0 (5.23) 

for any compact K C and any G G G))°(0, cxo), cf. P Formula (5.20)]. Moreover, 
seeing that $£ = Qg^QeUe) we may use estimate (15.11) to deduce from (I5.22p that 

Qe{.QeUe) 0 in Lf„^((0, T) x iF; R^). (5.24) 

Relations (I5.9p . (I5.24p imply the convergence claimed in (I2.2np . and, accordingly, 
(15.7p . In order to complete the proof of Theorem 12.21 it remains to show (15.Op . 
To this end, we will use the convergence result fl5.23p together with some uniform 
bounds derived in the following section. 

5.5 Uniform pressure estimates 

Our ultimate goal is to show fl5.6p . Obviously, using the uniform bounds fl3.14f[T25p . 
the convergence fl5.23p implies that 

^ ^ 0 in L^(0,T; L^gg^j^(JF)) for any compact K <Z VL. (5.25) 
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Consequently, in order to deduce fl5.6p we have to establish; 


(i) spatial compactness of 


^?£0e Q 
eg 


, . 1 -i- r 1 • ~ 

(ii) better mtegrabihty of the expression - - ---. 

As for the spatial compactness claimed in (i), we refer to Masmoudi [T^ Propo¬ 
sition 4.1], Specihcally, denoting {a; 5 } 5 >o the family of regularizing kernels, it can 
be shown that 


Us * 


^?e 0 £ 0 


Qe^e - Q 


LP{0,T-,L^{K)) 


—)■ 0 as h —)■ 0, 1 < p < cxo, (5.26) 


uniformly for £ > 0 for any compact K C hi, which, together with fl5.25p gives rise 
to 

^ 0 (strongly) in L^((0,T) x K)) for any compact K C hi. (5.27) 

As for (ii), we can proceed as in 0 Section 6 , formula ( 6 . 8 )] to show that 



{pe©e<p/2}niC 


(^£ 0 ^)'^+^ dx dt + 


'0 J{gsee>2g}nK 


(PeGe)'^’*'^ dx dt < e^c{K) 


(5.28) 


for any compact K <zVt. 

We can now combine (15.271) . (I5.28p to show (15.6p on compact subsets of hi. This 
can be shown by repetition of the argument from 0 Section 6 , formulas (6.5) and 
(6.9)] when is replaced by Qs^e- We have proved Theorem 12.21 


6 Concluding remarks 

The rather inconvenient restriction 7 > 3 in Theorem 12.21 is purely technical and 
has been effectively used only in Section 15.51 more specihcally, only in the proof 
of spatial compactness fl5.26p . It could have been relaxed should we have more 
information on integrability of An alternative approach to the problem would be 
to use the relative energy inequality, similarly to im. Such an approach, however, 
would require Strichartz type space-times estimates for the operator A in hi. To 
best of our knowledge, validity of these estimates is not known and maybe even not 
true in general, at least for the unbounded slab hi = x ( 0 , 1 ). 

Our hnal remark concerns existence of the weak solutions enjoying the regularity 
required in the present paper, in particular satisfying 02.41) or, alternatively, 02 . 6 p . As 
already mentioned in the introduction, an existence result in the variables [g, u, Z = 
p0] has been established in [1^. Next, we claim the following result that may be 
seen as a straightforward consequence of the theory developed by DiPerna and Lions 
0 . 
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Lemma 6.1 Suppose 


g > 0, Q E L^((0, T) X 12), Qg < Z < Qg, 0 < 0 < 1 < 0, 


and 


dtg + div 3 .(^)tt) = 0, dtZ + dS.Yx{Zu) = 0 m (0, T) x 
where the velocity field belongs to the class 

Then 


Zjg for g> 0, 


0 = 


G L°°(( 0 ,T) X n) 


1 for g = 0 
satisfies the transport equation 

dtQ + u- V^0 = 0. 


( 6 . 1 ) 

( 6 . 2 ) 

(6.3) 


(6.4) 


Proof: 

Using the regularizing procedure of DiPerna and Lions [ 8 ] we obtain 

dtps + u-Vxgs + gsdiVxU = r}, (6.5) 

dtZs + u ■ VxZs + ZsdiVxU = rj, ( 6 . 6 ) 

where, in view of fl 6 .ip . fl6.3p and Friedrich’s lemma, 

—)■ 0 in L^((0,T) x M^) as <5 —)■ 0. 


Next, we multiply fl6.5p by —(Z^ + A)/(p 5 + A)^, fl 6 . 6 l) by (p 5 + A)“^, respectively, 
where A > 0 is a positive constant. Summing up the resulting expressions we deduce, 
after a straightforward manipulation. 


dt 


f Zs + aa 

\ + A / 


+ div^, 


(£i±A). 

\ + A / 


{Zs + A)p5 ^ A 

(Pi + A)^ + A 


div3;it 


_ 1 + A ^2 1 

By virtue of fl 6 .ip . we may let hrst 5 —)■ 0 and then A —?■ 0 to recover the desired 
equation fl6.4p . 


Q.E.D. 

The solutions constructed in im enjoy the properties fl 6 .ip . fl 6 . 2 p . in particular. 
Lemma 16.11 yields fl2.6p at least for the values of the adiabatic exponent 7 > 2. In 
order to prove existence for the full range 7 > 3/2, one has to introduce the “artihcial 
pressure” proportional to AZ^ similarly to [T^ Chapter 3] and then perform the limit 
(5 —)■ 0 adapting the technique of Michalek j21) . 
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